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^ ! Abstract: Approximate symmetries of geodesic equations on 2-spheres are stud- 

OO '. ied. These are the symmetries of the perturbed geodesic equations which represent 

approximate path of a particle rather than exact path. After giving the exact sym- 
metries of the geodesic equations, two different approaches to study the approximate 



m 



in 

O I symmetries of the approximate geodesic equations show that no non-trivial approx- 

imate symmetry for these spaces exists. 
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1. Introduction 

For nonlinear problems, analytic solutions are rare and hard to obtain. Lie group 
theory provide a systematic and unified approach in search of analytical solutions 
[1, 2, 3]. Instead of solving the nonlinear equations directly, a system of over deter- 
mined linear equations is studied using the theory. Once the so-called symmetries 
of the equations are calculated, similarity solutions can be produced. By defining 
canonical coordinates, it is also possible to transform the equations into a much 
simpler form. Another approach allows to find approximate symmetries of approxi- 
mate differential equations. Generally, the perturbed term in a differential equation 
corresponds to some small error or correction. Therefore the resulting equations are 
defined approximately depending on a small parameter and such equations occur fre- 
quently in applications. The theory of approximate groups provides a regular method 
of calculating the perturbation directly, without using the complicated group trans- 
formations. The calculation is based on approximate symmetry equations. Recently, 
approximate symmetries have been used to find the approximate solutions of some 
partial differential equations [4, 5] as well. 

Let us consider an approximate (perturbed) equation in variables t and y, with 
a small parameter e, i.e. 

E{t,y,e)^0. (1.1) 

We write it as 

E{t,y,e) = Eo{t,y) + eE,{t,y) ^ 0, (1.2) 

where Eq and Ei are respectively the exact and approximate parts of the approximate 
differential equation (1.1). For (1.2) the infinitesimal generator can be written as 

[1,2] 

X = Xo + eXi, (1.3) 

where 

^ . d d 

at oy 

^ . 9 d 

Here ^o, G; Vo ^''^^ Vi ^^^ ^ill functions of t and y. Xq is the generator of the exact 
differential equation Eq, and Xi is the generator of the approximate differential 
equation Ei. Hence we can write 

X = (Co + e6) ^ + (^0 + erji) ^. (1.4) 



More generally 



where 



^ = 6 + 6^1, 

T] = r]o + er]i. 
A second prolongation [1, 2] is given by 

X = e^ + r/l- + r7« A + ^(2)^. (1.6) 

ot oy oy' oy" 

Equation (1.2) is approximately invariant under the approximate group of transfor- 
mations with the generator given in (1.3) if and only if Xi^l^-o = o(e) [2], or 

[Xoi^o + e(Xoi?i + Xii?o)]|^,o = «(e)- (1-7) 

This is the determining equation and X is called an infinitesimal approximate sym- 
metry or an approximate symmetry or an approximate operator [2], admitted by (1.2) 

If Xq is a generator of the unperturbed differential equation 

Eo = 0, (1.8) 

obtained by solving the determining equation 

Xo^o|^,o = 0' (1-9) 

then we define the auxiliary function H by 

i7 = ^Xo(i?o + ei?i)|^„^.,^^=o, (1.10) 

and find an approximate symmetry, (1.3), of the perturbed differential equation (1.2) 
by solving for Xi[2] in 

XiEo|£;«0 + ^ = 0- (1-11) 

Note that (1.11), unlike the determining equation (1.9) for exact symmetries, is 
non- homogeneous. 

The geodesic equations on a manifold can be regarded as a system of second 
order ordinary differential equations. The curvature of surfaces can be positive, zero 



or negative. Symmetries of geodesic equations for these three types of manifolds have 
been studied [6]. Interesting connection between these symmetries and the isometries 
[7] of these spaces have also been found. In this paper we explore approximate 
symmetries of geodesic equations of 2-spheres. Perturbations in geodesic equations 
correspond to approximate path of the particle as opposed to the exact one. It may 
be noted that we study approximate geodesies on a 2-manifold rather than exact 
geodesies on a perturbed manifold. 

2. Approximate symmetries of geodesic equations on a 2- 
sphere 

A geodesic is a locally length-minimizing curve and give the shortest distance between 
two points. Equivalently, it is a path that a particle which is not accelerating would 
follow. The geodesic equations on a manifold can be regarded as a system of second 
order ordinary differential equations. In coordinates x"" these can be expressed as 

[8,9] 

x" + Tlx'^x' = 0, (2.1) 

where dot represents derivative with respect to the arc length parameter s and F^^ 
denote the Christoffel symbols, which for a metric Qab, are defined by 

nc = 2^"' idbcd + 9M,c - 9cd,b) ■ (2.2) 

Here (/"^ denotes the inverse of the metric and the range of the indices gives the 
dimension of the space. 

Exact symmetries of the unit sphere have been given in Ref. [6]. We first 
outline the derivation of these symmetries and then use them for constructing the 
approximate symmetries. The metric for the sphere can be written as 

ds'^ = d9^ + sin^ 9d(j)'^ . (2.3) 

The geodesic equations (2.1) for this metric are given by 

El : ^'-sin^cos^^^ = 0, 

^2 : + 2cot^^0 = O. (2.4) 

These are second order ordinary differential equations. We apply the second prolon- 
gation [6], 

^^ ^ d , d o d 1(9 n(9 . d o d ,^^, 

^ds ' de ' d(f) '''de ' d(j) '89 '90 



of the symmetry generator given by (1.3), to both the geodesic equations. Here C,, 
1]^ and 77^ are all functions of s, 6 and 0; ri\ and rj'^g are functions of s, 9, 0, 6 and (p; 
and ri\g and rj'^gg are all functions of s, 6, 0, 6, 0, 6 and 0, Then 



Xi?! I Ei=0=E2 — 0) 
X-E'2 I Ei=0=E2 = 0) 



respectively yield [6] 



'q-^gg — 2 sin ^ cos 6(j)r]'^s — rj^ (cos^ 6 — sin^ 6) 0^ 
ryf,, + 2 cot ^ Ur]^^ + 9r]l) - 2^0 csc=^ 9r]^ 



\ei=0=E2 ~ ^' 



If we write 



then 



D 



d ■ d ■ d • ■ (9 

7- + 0—- + 0-— + 9^ + 0^ 
as 89 ^50 a^ ^90 



£;i=0=_B2 



0. 



vis = Dr^' - 9Di = nl^-9 {4 - i,) - 9^ie - 

r,l = Drf - 0De = vl + ^ {r,l - Q + 9ril - 4>% - ^06, 



vis 



Dv!s - m 

Vl + 9 {2r]ls - Css) + 9' {rjl, - 2^,o) - 9^iee + 24>iil^ 
+02ryJ^ + 2^0 (77!^ - is^) - 29^ie4. - 0^%4> 
+9 (77,1 - 2e. - 3^6) - 2^0e<^ + 4>r]l - 9k<p, 



V%s 



Dvl - (f>m 



2(^2 
2, 



+ OVo 



= vl + 29r,l + 9W,s + (2^'^ - U) + 
-^^U^ + 2^0 {r,l^ - ise) - 9'kee - 

-Oke +'<p(ril- 26 - 30^^) - 2^06- 
Inserting (2.9) - (2.12) into (2.6) and (2.7) yields 

[nl + 6 {2rile - iss) + 9^ [vie " 2^^) - O'iee + 2077^^ + 0^77^^ 
lU - U) - '^^''ke<t> - ^'/''U + ^ [vl - 26 - 3^6) 
+ H - (^k<p\ - sin 2^0[77,2 + (775 - 6) + er,l - ^U 

- r]^ cos 29^1 Ei=o=E2 = 0, 



(2.6) 

(2.7) 

(2.8) 



(2.9) 
(2.10) 



(2.11) 



(2.12) 



(2.13) 



m' ■■ vis 

vis = 0, 


= 0, 


■■ He - iss 
vie + cot 9r]l 


= 0, 
= 0, 



[Vl + 29r^% + eWee + 4> {H^ - Q + 0' {vl^ - 2^^^) - ^ 
e<p - Cse) - 0'kee - 29<P%^ + Oril - O^e + 4> {vl - 2^. - 30^^] 

-2eke] + 2 cot 9i<PH + d [vl " Q - O'Ce + <Pvl " OM 
+Hvl + Hvl- e.) + Ov'o - 4>% - dke}) - 2^0 csc^ ev%,=o=E, = 0. (2.14) 

Inserting the values of 6 and from the geodesic equations (2.4), and then comparing 
the coefficients of the powers of 9 and <p, we obtain the following system of partial 
differential equations. 

2.15) 
2.16) 



2.17) 
2.18) 

2.19) 
2.20) 

2.21) 

2.22) 

2.23) 
2.24) 

2.25) 
2.26) 

2.27) 
2.28) 

2.29) 
2.30) 



<P:2vl^-sm2erj^,=0, 
2l^%-^ss + 2cot9r]l = 0, 

9^ : vis - 2U = 0, 
rils + 2coi9ril = Q, 

h"^ : ri^^ + sin 9 cos 9rjl — rj^ cos 29 — sin 29r]'^ = 0, 
-^ 2^^,^ + sin 9 cos 9r]'^ + 2 cot 9r]j, = 0, 

9^ : Cee = 0, 

- sin ^ cos ^(^61 = 0, 



■ 2 ivl^ - 60) - sin 29r]'s - 2 cot 9^1 = 0, 
{Ve^ - Cse) + cot 97]l - csc^ 9r]^ = 0, 

9^ '■ —i<t><t> ~ sin 9 cos 9^0 = 0, 
-2^e^ + 2cot^e0 = O. 



Equations (2.25), (2.26) and (2.29) yield 

e = bi(s). (2.31) 

Substituting the value of ^ into (2.21) and solving with (2.15), we get 

V' = [fMs + /2(0)] + /3(0)s + fM- (2.32) 

From (2.16), (2.18) and (2.22) we obtain 

v' = -coteM) + M^). (2.33) 

Substituting the values of rj^ and r/^ into (2.19) we get 

V' = [As + hm e + Bs + U{<P). (2.34) 

Then (2.17) yields 

^ = As^ + ciS + Co. (2.35) 

Now (2.23) and (2.27) give 

e = cis + co, (2.36) 

r]^ = C3 cos + C4 sin (/), (2.37) 

rf = cot 6 (c4 cos — C3 sin (/)) + C2- (2.38) 

Therefore, the exact symmetries come out to be 

Y-^Y- ^ ^ - ^ 

-^0 — 7^1 ^1 — "StT'^S — 77-5 
OS OS 0(p 

d d 

X3 = cos (p— - cot e sm 0—, 

d d 

X.4 = sin (f)—- + cot 9 cos (p—- . (2.39) 

00 0(j) 

In the next sections we discuss two different approaches for finding approximate 
symmetries of this 2-manifold. 

2.1 Approximate symmetries: First approach 

Here we will convert the geodesic equations (2.4) of the sphere into perturbed equa- 
tions by adding a general function and then try to find the approximate symmetries 
of these equations. 



The second prolongations of the five exact symmetries (2.39) of the geodesic 
equations of the sphere are given below. 

(i) XO = A 

Here ^ = 1, t]^ = = rj'^, therefore, equations (2.9) -(2.12) give ij^g = = 'q'^g, 
v\s = = rj'^gg so that the second prolongation becomes 

XO = ^. (2.40) 

(ii) X^ = .f 

Here ^ = s, r]^ = = r]'^, therefore, equations (2.9) -(2.12) give ij^ = —6, 
rfg = —0, ri\g = —29, rfg^ = —2(f), so that the second prolongation becomes 

x' = sl-el-^^-2el-24>t. (2.41) 

ds de d<j) de d(p 

(iii) X^ = & 

Here ^ = 0, t]^ = 0, t]"^ = 1, therefore, equations (2.9) -(2.12) give ^^^ = = ?]^^, 

r]\g = = 'q'^gg so that second prolongation becomes 

X= = |. (2.42) 

(iv) X"^ = cos (p-^ — cot 9 sin 0^0. 

Here ^ = 0, 77^ = cos0, 77^ = — cot 6* sin 0, therefore, equations (2.9) -(2.12) give 

vi = -0sin0, 



Tj'^g = 9 (sin (p csc^ ^) + (— cot 9 cos 0) , 



rj\g = —<p cos — sin 0, 

7^2^^ = 9^ sin (-2 csc^ 9 cot 6i) + 0^ cot 6^ sin + 26^0 (cos csc^ 9) 
+9 csc^ ^ sin + (— cot 9 cos 0) , 

so that the second prolongation becomes 

d ' / ' ' \ 

X^ = COS 07" cot 9 sin 0— sin — r + ( ^ sin csc^ 9 — (p cot ^ cos ) — - 

89 dcp d9 ^ ) (90 

— ( (p? cos + sin j ^ + [—29'^ sin csc^ 9 cot 9 + (j? cot 6* sin 

+26*0 cos csc^ 6* + 6* csc^ 6* sin - 0cot 6'cos0)^t. (2.43) 



(v) X^ = sm(f)i, + cot9cos(l)4i 



Here ^ = 0, 77^ = siiKJ), rf = cot6'cos(/), therefore, equations (2.9) -(2.12) give 
rj^g = (pcoscj), 

r}\ = —6 csc^ 9 cos 4> — (p cot 6 sin </), 
rj^gg = —0^ sin (/) + cos cp, 

7^2^^ = ^2 ^2 csc^ 6* cot 9 cos 0) + 0^ (- cot ^ cos (p) + 2^0 (csc^ 9 sin (/>) 
— ^ csc^ 9 cos (f) — (j) cot ^ sin <;/>, 

so that the second prolongation becomes 

d 3 • d / ' ' \ 3 

X^ = sin (p— — h cot 9 cos (p— — \- (p cos (p — r — ( 9 csc^ ^ cos + cot ^ sin ) — - 
39 d(p 39 \ J 3(f) 

/ •■ • \ 3 • • 2 

+ ( COS — 0^ sin j ^ + (2^^ csc^ ^ cot 9 cos H — cot ^ cos 

+2610 csc^ 61 sin - 61 csc^ ^ cos - cot 9 sin 0) ^r . (2.44) 

50 

Now, we approximate the geodesic equations of a sphere by using two arbitrary 
functions f{9), g{(p) and a small parameter e as 

El: 9- sin ^ cos ^0^ + e/(6') = 0, 

^2 : + 2cot^^0 + e^(0) = 0, (2.45) 

where 

El : 6'' - sin ^ cos ^02 = 0, 

£;^ : + 2cot^^0 = O, (2.46) 

are the exact geodesic equations. For Case (i) the auxiliary functions become 

Hi = ;X°(^--sin^cos^02 + e/(^))|^_,_,^.^,^(,)^, 
I f 3\ ,, . „ -2 



- (0) = 0, 



H2 = -X° (0 + 2cot^^0+e^(0)) |^+2eotee0+e5W=O 

" 7 (ai) (^ + 2cot^^0 + e^(0)) |^+2cotee^+.5W=o 

= i (0) = 0. (2.47) 



Hence Hi = = H2, so we cannot proceed further. 

Now we take Case (ii). Here the auxihary functions become 

Hi = -X^(^'-sin^cos^(/)2 + e/(^))|^_,i,,,,„,,^.+,^(,)^o 

1/(9 -d ■ d •■ a yd 
e\ds dd (90 86 d(j) 

X (^9 - sin ^ cos ^(^2 ^ gj(^)j le_sin0cose02+,/(e)=o 

= i (-2^- + 2sin^cose(/)2 + e(0)) b_,i„,eose<^^+./W=o 



H2 = -Xl (^0 + 2 cot 66^ + €g{<p) 



l0+2cote0</>+eg((^)=O 

e\ds dO 84) dO dcp 



X i^(t) + 2 cot 99(1) + eg{(P)j |^+2cotee<^+.9{0)=o 
= - (-24) - 2 cot 994 - 2 cot ^^0 + e (0) 



-^{-2{-eg{4))) = 2g{4). (2.48) 



Solving for X\ in 



Xi£^0 I £;l=o + -f^i = 0, 

XiE^ I „2=o + ^2 = 0, (2.49) 



gives the same equations as (2.17)-(2.30) with a change in (2.15) and (2.16) given by 

Vi + 2f{9) = 0, (2.50) 

7^^ + 2(7(0) = 0. (2.51) 

Equations (2.25), (2.26) and (2.29) yield 

e = 6i(s). (2.52) 

Substituting the value of ^ into (2.21) and solving with (2.50), we get 

r^^ = f,[9,4^)s + U9,4)-f{e)s\ (2.53) 



10 



Integrating (2.53) twice with respect to 6, and using (2.21), we get 

rile = fieeiO, 0)s + f2ee{0, 0) - fe9{0)s' = 0. (2.54) 

Comparing the coefficients of the powers of s, we obtain 

fwe = 0,f2ee = 0,fee = 0, (2.55) 

which imphes that /i, /2 and / are hnear in 9. Hence we can write 

f(9) = ce + d, (2.56) 

so that (2.53) becomes 

v' = [fii<P)s + /2(0)] e + M<j))s + fM) - {ce + d) s\ (2.57) 

From (2.51), (2.18) and (2.22), we obtain 

rfs = 0, (2.58) 

therefore (2.51) gives 

Vi + 2g{cp) = 0, (2.59) 

so that 

g{<P) = 0, (2.60) 

r/' = -cot^/6(0) + /7(0). (2.61) 

Substituting the values of t]^ and r/^ into (2.19) we get 

vU = 0,fii<P) = A,M^) = B, (2.62) 

which gives 

7]^ = [As + f2{(p)] e + Bs + U{(j)) - {ce + d) s^. (2.63) 

Thus (2.17) yields 

2 
bi (s) = As^ - -cs^ + as + Co, (2.64) 

so that (2.52) becomes 

2 
C = As^ cs^ + CiS + cq. (2.65) 
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Now, inserting the values of 77^ and rf in (2.23) and (2.27) gives 

/^(0) - cot ofM - cot efi{^)e - cot ^/^(0) = 0, 

f2{4>)d + /r(0) + ^sAsin2^ + i sin 2^/2(0) - ]^s\sin2e 

-s9A cos 29-9 cos 26'/2 (0) - sS cos 2^ - cos 2^/4(0) 

+ cos 29 {c9 + d) s^ + 2 cos^ ^/^(0) - sin 2^/^(0) = 0. 

Comparing the coefficients of the powers of 9 and s, we obtain 

c = 0,d = 0,/(^) =0. (2.66) 

From equations (2.60) and (2.66) it is clear that, for this case, we do not have any 
new symmetry. Similarly, we see that Cases (iii), (iv) and (v) also do not give any 
non-trivial symmetry. 

2.2 Approximate symmetries: Second approach 

Now we adopt another approach and take a more general function to make the exact 
geodesic equations (2.46) perturbed, and investigate the existence of approximate 
symmetries. Here we approximate the geodesic equations of a sphere as 

El : 9- sin 9 cos 9(j)'^ + e/(^, 0, 9, 0) = 0, 

E^ : ^ + 2 cot 99(j) + eg{9, 0, 9, 0) = 0. (2.67) 

As these are second order ordinary differential equations, we apply the second pro- 
longation 

,^^(9 ,8 r, d , d r, d . d r, d ,, 

^ds '86 ' d(P '''89 '50 '89 ' 8(j) 

where 

e = eo + eei, (2.69) 

r] = r]o + er]i, (2.70) 

and the infinitesimal generators, as before, are given by 

X = Xo + eXi. (2.71) 

Equations (2.67) are approximately invariant under the approximate group of trans- 
formations with the generator given in (2.68) if and only if 

XE" I E^=0=E'2 = 0, 

XE^ I ^^^^^^, _ 0, (2.72) 
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e-smecose<j)'^+ef{e,(j>,e,,j})=o 



0, 



4>+2cotee,j}+€g{e,(f),e,,j})=o 



0. 



so that we have 

:S.(e- sin e cos e^^ + e/(^, 0, ^, 0) 

X (^0 + 2 cot ee^ + eg{e, 0, ^, 0) 
Now, let us take the following form of the functions / and g introduced in (2.67) 
f{9, 0, 9, 0) = A:i + k29 + h(f) + h9 + h^ + ^9^ + A;70', 

g{9, 0, ^, 0) = /li + /l2^ + /l30 + hi9 + /l50 + /l6^^ + /i70^ 

where fci, ..., k-j, /ii, ..., /17 are all constants. Then (2.73) yields 

[r]\^ — 2 sin ^ cos 9(j)rj\ — rj^ (cos^ 9 — sin^ 6') 0^ + e{k2r]^ + A;3?7 



(2.73) 



(2.74) 



[77 



2 

ss 



+k^ri]^ + A;5?7' + 2k^9ri]^ + 2A;70r7y]|^^^o^^^ = 0, 

2 cot ^ ( 0?7i + ^r]^ ) - 2^0 csc^ 9r]^ + e(/i2?7^ + h^rf 

0. 



+/i4r7i + h^nl + 2/i6^r7i + 2/i70?7i)] U^=o=e. 
Inserting (2.9) -(2.12) into (2.75), we obtain 

[ril + 9 (277I, - U) + ^' He - 26e) - ^'6e + '^Ws^ + 0'^J<^ 
+2^0 [r,]^ - e,^) - 2^20^^^ _ ^2^^^ + ^ (^^1 _ 2^^ _ Wi,^ 

+ hl- OkA - sin 29^[ril + (77^ - Q + ^r^,^ - 02^,^ 



(2.75) 



r]^ cos294P'\e^=o=E2 = 



^2 ^-2 . 2es<^) - 0'e00 + 2^0 (?7e0 - 6e) 
72-2e.-30e<^)-2^0e,]+2cot^ 



[vl + 2^77^, + 9'r]ls + (277,2^ - Q + 

iHvl + {vl - ^s) - 9% + H - HU} + oWs +Hvl- Q 

+OVe - 4>% - dke}) - 2^0 csc^ Or]\^^,^^^ =(0.76) 

Inserting the values of 9 and from the geodesic equations (2.67), the above equations 
take the form 

l^i + 9 (277I, - Q + 9' [vie - 26e) - 9%e + 2077^^ + 0^77^,^ 

+2^0 (?7,^^ - ^s<t>) - 29^Ce<^ - 9kiH + (sin 9 cos 9^ - tf) 

{ril - 2is - me) - 2 (sin 9 cos 9^ - ef)k<t> - (2 cot 99^ + eg)r]l 

+9{2 cot 99^ + eg)i^ - sin 2^0[772 + (77^ - ^ + 9^ - ^U " ^e 

-ri' cos 29k + <hv' + hv^ + {h + 2k^9){r,l + 9 {r,l - Q - 9^ie + 0^7^ 

+ (^5 + 2A;70)(r72 + (775 - 6) + %' - kU ' ^e) = 0, 
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I2 



{Ve<t> - ise) - 0'^(j)iee - '2d9Ce<p + (sin 9 cos 6*0^ - ef)i]Q - (sin 9 cos i 
-ef)ke - (2 cot 99<j) + e^) (vl - 2^, - 30C<^) + 29(2 cot ^^0 + e^)^^ + 2 cot ^ 

iHvl + [4 - Q - e'ie + hi - HU) + e{vl + H^l- Q + (^4 

-<i>% - Oke]) - 29<i)csc^ 9r]' + e{h2v' + hsV^ + {h^ + 2h(,e){ril + 9 {ql - Q 
-d% + kl - Ok^) + (^5 + '^hj^Ws + (^J - 6) + 04 - 4>% - Oke) = C2.77) 

Now, we use (2.69) and (2.70) in (2.77), and separate terms with e° and e^, and 
ignoring terms with e^. For e° we get 

VL + {2vle - Coss) + 9' {v'oee ' 26e) - O^Coee + 2(/>r/o% + kvl^^ 

+29^ [rjlo^ - ^os^) - '29^koe<p - OkCo4>4> + sin 9 cos 9kiVoe 

-2^s - ^koe) - 2 sin ^ cos 9k^ocf> - 2 cot 99^1^ + 2 cot 99^ko<i> 

- sin 29^[vl + 'P (^00 - Cos) + He " kio<t> - Okw] - vl cos 29^ = 0, 

vL + 2^%2^e + OWoee + (2^'<^ - Coss) + </>' (??o^^ - 2^,.^) - kCo^^ 
+29<j) [rjQg^ - ^ose) - 9^koe9 - 29kiw<t> + sin 9 cos 9krile - sin ^ cos 9kiw 

-2 cot ^^0 (r/2^ - 2e. - 30eo<^) + 4 cot ^^^^^^^ 

+2 cot 9{kvls + ^ (%'e - eos) - ^'eoe + 0%'^ - 9ko^] 

+d{rils + {vl^ - eos) + ^%'e - kio4> - Okw}) - 2^0csc2 9r,^ = 0.(2.78) 

Comparing the coefficients of the powers of 9 and cp, we obtain a system of partial 
differential equations, which on solving yields the following solution 

^0 = cis + Co, (2.79) 

r^Q = C3 cos + C4 sin (/), (2.80) 

ril = cot 9 (c4 cos (/> — C3 sin cp) + 02- (2-81) 

We follow the same procedure for e^, compare the coefficients of the powers of 9 and 
(p, and obtain a system of partial differential equations. We denote 

Cos = A, Til = B, rjl = C, r,l^ = D, r,l^ = E, r^l, = F, (2.82) 



14 



so that the system of equations become 

{e^f : ri\^^ + 2{h + he + h(j))A - (/ii + h^d + /i30)D + hB + hC = 0, 
r^^^ - {h + A;2^ + A;30)F - (/ii + /i2^ + h^){E - 2 A) + h^B + h^C = 0, 

^ : 2r]le - ^iss + hA - KD + A:5F = 0, 
2r]l^Q + 2 cot er]l^ - hF - h^E + h^A + /15F = 0, 

: 277i% - sin 29ril^ + A;5A - h^D + A;4Z^ + hE = 0, 
2??L<^ - ^iss + 2 cot 6'?7i\ - hF + hA + hD = 0, 

^' : vlee " 2^1^^? " /^6^ = 0, 
?7i0e + 2 cot 9r]lo - hF - hE = 0, 

0^ : r]l^^ + sin 6* cos 9r]lg - ri\ cos 26* - sin 20iril^ - /i7£) + 2hE = 0, 
^100 - 26s0 + sin ^ cos 9r]fe + 2 cot 9r]l^ - /cyF + /ly^; = 0, 

9 : Gee = 0, 
^^ : -600 - sin6'cos6'Ge = 0, 

^0 : 2 (?7|,^ - eis^) - sin 29r]ls - 2 cot 9r]l^ + 2hD + 2hF = 0, 
2 (^?e0 - ^ue) + 2 cot 9r]lo - 2 csc^ ^77^ + 2hD + 2/17^ = 0, 

04>^ : -6</.</. - sin 6' cos ^^0 = 0, 
-2^w^ + 2 cot 9^1^ = 0. 



2.83) 

2.84) 

2.85) 
2.86) 

2.87) 
2.88) 

2.89) 
2.90) 

2.91) 
2.92) 

2.93) 
2.94) 

2.95) 
2.96) 

2.97) 
2.98) 



Equations (2.93), (2.94) and (2.97) yield 

Ci = bi{s). (2.99) 

Substituting the value of C, into (2.89) and solving with (2.83), we get 

h = 0, (2.100) 
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_1a;25s2 + ((^i(0)^ + g,{cP))s + (73(0)^ + <74(0) + l^D^ , (2.101) 

where gi{4>)-, 92(4')^ 9z{4>) ^"^^ 9i{4>) ^-^e functions of integration. From (2.84) and 
(2.86), we obtain 

ki = k2 = h = hi = h2 = h = 0, (2.102) 

Vf = 2^[- cot e{hF + hiE - hiA - h^F) - hFe - h^Eg + h,,Fe\ + 73(6', 0), (2.103) 

where /3(^, 0) is a function of integration. Using the above value of rjl in (2.90), we 
get 

k4 = h4 = h5 = 0, (2.104) 

r^2 = (6'cot^-lnsin^)A;6D- -6'/i6l^ - cot ^ 6/5(0) + t/6(0), (2.105) 

where 5'5(0) and (76(0) are functions of integration. Using the values of 'ql and rjf in 
(2.87), gives 

9i{<P) = C5, 

^2(0) = --A;5(C + A0) + C6. 

Now (2.88) implies 



Similarly (2.95) and (2.91) give 



h = o, 


(2.106) 


C5 = 0, 


(2.107) 


^ = C7S + Cs- 


(2.108) 


h = o, 


(2.109) 


5^3(0) = A;6-B + Cg, 


(2.110) 


A;6 = 0, 


(2.111) 


cg = 0. 


(2.112) 


Cii sin (A;7 tan 6F + /17-B) , 


(2.113) 



kr 
5'5(0) = cio sin — cii cos — — tan 9 csc^ ^D 

-^S - ^0D + fc7tan^F(^ + 1), (2.114) 
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fi'400 — COS 2^(74 + 2 cos^ 9g^^ — sin 2^*76^ — h^D + 2kjE = 0. (2.115) 

Again, using the values of 'ql and rif in (2.96) and (2.92), we get 

96{(p) = ci2, (2.116) 

/i7 = 0, (2.117) 

kr = 0. (2.118) 

Therefore, we have 

6 = C7S + C8, (2.119) 

^1 = cio COS0 + cii sin0, (2.120) 

tjI = cot 9 (cii cos (p — cio sin 0) + c, (2.121) 

which are trivial, as all the constants ki, ..., kj, hi, ..., hj in (2.74) become zero. Hence 
we see that this case also does not give any non-trivial approximate symmetries. 

3. Conclusion 

One of the methods to solve differential equations is the symmetry method by which 
we can solve differential equations, reduce the order of an ordinary differential equa- 
tion and can reduce the number of independent variables in a partial differential 
equation if it is invariant under a one-parameter Lie group of point transformations. 
The method of 'canonical variables' is one of the basic procedures for the integration 
of ordinary differential equations with known symmetries. Conservation laws for the 
Euler-Lagrange equations can be constructed when their symmetries are known and 
hence Lagrangians can be constructed. 

For tackling differential equations with a small parameter, the method of 'approx- 
imate symmetries' can be used. Considering the geodesic equations of a sphere as a 
system of two ordinary differential equations we have investigated their approximate 
symmetries. Two different approaches have been adopted to find these approximate 
symmetries. Firstly, by converting the geodesic equations into perturbed equations 
by adding the general functions f{6) and g{<f)) and taking a symmetry of the geodesic 
equations of the sphere. In the second approach we use the functions 

f{9, 0, 9, 0) = A:i + k29 + A;30 + kj + k^^ + k^O^ + k^^'', 
g{9, 0, 9, (j)) = hi + h29 + /i30 + hi9 + h^^ + /ig^^ + /l70^ 

to convert the geodesic equations into perturbed equations. In both the cases, no 
non-trivial approximate symmetry has been found. 
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It would be interesting to see if one can obtain non-trivial approximate symme- 
tries for these manifolds by some other method. 
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